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Differential evolution (DE) algorithm has been applied as a powerful tool to 
find optimum switching angles for selective harmonic elimination pulse 
width modulation (SHEPWM) inverters. However, the DE’s performace is 
very dependent on its control parameters. Conventional DE generally uses 
either trial and error mechanism or tuning technique to determine appropriate 
values of the control paramaters. The disadvantage of this process is that it is 
very time comsuming. In this paper, an adaptive control parameter is 
proposed in order to speed up the DE algorithm in optimizing SHEPWM 
switching angles precisely. The proposed adaptive control parameter is 
proven to enhance the convergence process of the DE algorithm without 
requiring initial guesses. The results for both negative and positive 
modulation index (M) also indicate that the proposed adaptive DE is superior 
to the conventional DE in generating SHEPWM switching patterns. 
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1. INTRODUCTION 

Selective harmonic elimination pulse width modulation (SHEPWM) is one of the methods to 
generate sinusoidal voltage output from single- or three-phase inverters. However, the SHEPWM often 
experiences with the complexity of calculation the switching angles which take quite long computing time to 
achieve the optimum values. That will be the disadvantage for on-line and real-time applications [ 1 ]-[3]. 

Evolutionary algorithms such as genetic algorithm (GA) [4]-[6] and differential evolution (DE) 
[7]-[ 10] have been used to reduce the complexity of the SHEPWM problems which initial values are usually 
not guessed properly. The evolutionary algorithms can generate a large jump step to find the most 
approximate optimal solutions. However, the runtime behavior of evolutionary algorithms associate with 
control parameters, i.e., crossover (CP) and mutation factor ( F ). The control parameters are very sensitive. 
Choosing the proper parameters is quite difficult and dependent on the problems [11], Choosing 
inappropriate value of CR and F will make the search process suffer from slow convergence speed, 
premature convergence and stagnation problems [12], [13], 

The objective of this paper is to propose an adaptive control parameter for DE in order to avoid or 
reduce requirements for tuning consumption time due to uncertainty of the CR and F values. As the adaptive 
function, a triangular probability distribution is employed to adjust the control operators of DE algorithm. 
The solution patterns for both positive and negative values of modulation index ( M ) in the range of the 
SHEPWM capability are also investigated. Additionally, the performance of the proposed DE is compared to 
the conventional DE. 


Journal homepage: http://iaesjournal.com/online/index.php/UPEDS 








1382 □ 


ISSN: 2088-8694 


2. SELECTIVE HARMONIC ELIMINATION IN THREE-PHASE INVERTER 

Figure 1 illustrates a generalized PWM pattern. There are only odd harmonic components exist due 
to the symmetries in the PWM waveform. The Fourier series for this output voltage waveform is defined as 
follows: 

00 

v(cot) = sin(/7<rtf) (1) 

n=l,5,7,... 

The Fourier coefficients ( V„) are then given by 

N 

l-2^(-l) ; cos focij) 

H 





Figure 1. Generalized PWM waveform 


Equation 2 can be solved in order to obtain N switching angles, i.e., ai to <x v by equating N -1 harmonics to 
zero. In three-phase inverters the modulation index is set as M = VjIV^c and the normalized magnitudes of 
harmonics need to be eliminated are as follows [14]: 



All triple harmonics are absent from Equation 3, e.g., for N =9, there are 8 harmonics to be eliminated, i.e,/ 5 , 
/ 7 >/ii>/i 3 >/i 7 >/i 9 >/ 23 > an d fis- The objective function of the problem is 

f(a)=£ l +£ 5 +£ 1 +... + £ ON _ 2) (4) 
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3. ADAPTIVE DIFFERENTIAL EVOLUTION ALGORITHM 

Differential evolution algorithm utilizies control parameters, i.e., mution factor (F) and crossover 
rate ( CR ) to control perturbance and improve convergence of the optimization process. Both F and CR have 
values in the range of [0, 1]. The flowchart of proposed adaptive DE to optimize N switching angles is 
presented as in Figure 2. The first step of the flowchart is generating the first generation of population. 


a, 


(i) 


= ry ■ ■ 4 - 
^minj 1 


randj (a 


(5) 


where a mmj and a maxj are lower and upper bounds of the swithing angles, randj e [0, l],y = 1 , 2, i= 1 , 

2, NP, and NP is the population size. Evaluating the population is conducted to determine which 
switching angles of the current generation will be the best candidates in order to satisfy the constraints of the 
objective function. 


fbeJ" =/(«teJ 1> )> a, 


(i) 


bestj 


: a: 


(i) 


(6) 


New individuals are evolved by using mutation and crossover operators. The operators are employed to 
increase the population diversity and promote faster convergence. Mutation operation has a responsibility to 
create mutant individuals, while crossover operation is to create trial individuals. The individuals in this case 
are the switching angles. Both mutation and crossover processes are shown as in Figure 3. In the mutation 
process, the indices ra, rb, rc, rd, and re are five mutually distinct integers taken randomly from {1, 2, 3,..., 
NP}. Further, the integers i, ra, rb, rc, rd, and re must be different. In selection process as seen in Figure 4, 
the trial individuals (t,j) become the best next generation individuals if they offer equal or lower values of 
the objective function than those of their parent. 



Figure 2. The flowchart of adaptive DE 
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Figure 3. Mutation and crossover operations 



Figure 4. Selection operation 


3.1. Adaptations of Mutation and Crossover 

The adaptations of mutation and crossover are employed to avoid ineffective either trial and error or 
tuning processes in order to meet the best values of F and CR. The control paramaters at the Gth-generation, 
i.e., F, <a> and CR,' 0 ’ are adapted by using triangular probability distribution as follows: 


Z <G) I x m,n + ^ ran dj(X niax - X m J(X moi - X m J 


if randj < r 


X max ~4(!- randj )(X max - X min )(X max - Z mod ) otherwise 


'T ( X mod X min )/(X max -X min ) 


(7) 


The triangular probability distribution on [0, 1] is defined by three values, i.e., lower limit X mi „, 
upper limit X max , and mode X mod , where X min < X max and X min < X mod < X max . There are three possible cases [15] 
as shown in Figure 5, i.e., (i) right triangular (X mad = X min ), (ii) middle triangular (X min < X mod < X max ), and (iii) 
left triangular (X mod = X max ). 
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Figure 5. 



Tringular probability distributions 


4. RESULTS AND ANALYSIS 

The adaptive DE is applied to minimize the objective function. The proposed DE is implemented by 
using MATLAB software package in order to generate SHEPWM switching patterns. There are three cases 
of triangular (as in Figure 5), which are investigated to adapt both F and CR. 

4.1. Triangular Adaptive DE 

There are 9 possible combinations for F and CR. i.e., (F = right triangular, CR = right triangular), 
(F = right triangular, CR = middle triangular), (F = right triangular, CR = left triangular), ( F = middle 
triangular, CR = right triangular), (F = middle triangular, CR = middle triangular), (F = middle triangular, 
CR = left triangular), (F = left triangular, CR = right triangular), (F = left triangular, CR = middle triangular), 
and (F = left triangular, CR = left triangular). From the aforementioned combinations, there is only one 
combination is applicable to determine SHEPWM switching patterns, that is (F = left triangular, CR = right 
triangular). The other combinations will result in the non-corvergence DE. 

Figure 6 shows F and CR tracking of the triangular adaptive DE in order to satisfy the best 
switching angles for the SHEPWM with N = 9, M = 0.05, and f(a) < 0.0001. It needs 25 generations to meet 
the optimum switching angles as shown in Figure 7. The triangular adaptive DE is faster than the 
conventional DE using fixed F and CR due to it has a large jump to the optimum solutions. According to the 
report in [14], the best fixed F and CR values are 0.26 and 1.00, respectively. However, it needs 45 
generations to solve this switching angles problem, as shown in Figure 8. The final best switching angles to 
be found by both DE are a, = 11.7423°, a 2 = 12.0905°, a 3 = 23.7342°, a 4 = 24.1551°, a 5 = 35.7282°, 
a 6 = 36.2035°, a 7 = 47.7291°, a 8 = 48.2380°, and a 9 = 59.7398°. 
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Figure 6. F and CR adaptations 
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Figure 7. Triangular adaptive DE tracking the best switching angles 



Figure 8. Conventional DE tracking the best switching angles 


4.2. SHEPWM Switching Patterns 

There are 4 types of possible switching patterns for N = 9 with both negative and positive values of 
M, as presented in Figure 9. The upper panels correspond with the type-1 of both negative and positive 
switching patterns, and other panels correspond with type-2, type-3 and type-4, respectively. The negative 
and positive values of M are ranged between -1.5 < M < 0.0 and 0.0 <M< 1.5, respectively. The computation 
process is employed with a step of 0.05. The number of generations to generate every type of patterns is 
summarized as in Table 1. The number of generations is the averaged values over 100 runs and with 100% 
success rate. 
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Figure 9. SHEPWM swithing patterns for N =9 with M negative and positive 


As shown in Table 1, the adaptive DE is faster than the conventional DE in generating switching 
patterns of SHEPWM. Note that the conventional DE used fixed F and CR of 0.26 and 1.00, respectively. It 
is discovered that the triangular adaptive F and CR improve the speed of DE to reach the optimum values of 
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switching angles of the SHEPWM problems. The adaptive DE requires less number of generations when 
compared to the conventional DE in generating SHEPWM switching patterns. 


Table 1. Number of Generations of DE 




Number of generations 


Patterns 

Adaptive 

Conventional 



DE 

DE 


Type-1 

44 

82 

M 

Type-2 

45 

77 

negative 

Type-3 

47 

74 


Type-4 

45 

76 


Type-1 

91 

182 

M 

Type-2 

91 

197 

positive 

Type-3 

97 

160 


Type-4 

49 

77 


5. CONCLUSION 

The application of triangular adaptive control parameter has an advantage in increasing the DE’s 
performance to solve the SHEPWM switching angles problem. It reduces half of the number of generations 
required to meet the optimum values of switching angles in comparison to the conventional DE. It works 
relative fast with 100% success rate for both negative and positif values of modulation index. These results 
provide usefull information to help designers to make a decision in their SHEPWM inverter for real-time 
applications 
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